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The extended Dynamical Mean Field Equations (EDMFT) are analyzed 

using semiclassical methods for a model describing an interacting fermi-bose 
system. We compare the semiclassical approach with the exact QMC (Quan- 
tum Montecarlo) method. We found the transition to an ordered state to be 
of the first order for any dimension below four. 



The dynamical mean field theory (DMFT) , a recently developed many body approach to 
strongly correlated electron systems has been very successful in unraveling non perturbative 
problems such as the Mott metal to insulator transition.0 In spite of its many successes, 
this technique has several limitations resulting from its single site character and from the 
lack of feedback of the non local collective excitations on the one particle spectra. Several 
approaches are being pursued to extend the scope of the DMFT method. In this paper 



energy while incorporating feedback effects of the charge and spin dynamics in the one 
electron properties. This method gives rise to quantum impurity problems with fermionic 
and bosonic baths that need to be solved self consistently. This method has already been 
applied to wide class of models, such as the spin fermion model,i fermions interacting with 
long range (Coulomb) electron-electron interaction,! electron-phonon systems^ and frus- 
trated magnets.! 

The EDMFT equations are more involved than the conventional DMFT equations be- 



I. INTRODUCTION 



we explore an extension of the DMFT method 




cause they involve a solution of a self consistency problem in an additional bosonic sector 
and only recently a full numerical analysis of the self consistency conditions of EDMFT 
was carried out! The interpretation of the EDMFT instabilities is also not as straightfor- 
ward as in DMFT because bosonic and fermionic propagators involve very different regions 
of momentum space, and a formulation of EDMFT for ordered phases was only obtained 
recent lyi 

The purpose of this paper is to develop the EDMFT approach further by analyzing several 
aspects of this method: a) We implement a semiclassical technique for its solutionBffl 
and compare its results to the earlier QMC studyB to test its accuracy We show that 
the analytic treatment is in satisfactory agreement with exact (QMC) results in the high 
temperature regime of the three dimensional (3d) model and provides analytic expressions 
for various physical quantities, b) We extend this study to the case of two dimensional (2d) 
phonons, which had not been treated in ref.i We demonstrate that in the 2d case the EDMFT 
treatment at finite temperatures, if it produces an ordering transition it is necessarily of the 
first order. This analysis applies to a very general class of models including those used in 
ref.i 

We also analyze the EDMFT equations in the ordered phasejl for a simple spin model. 
This analysis clarifies the strengths and the limitations of the EDMFT approach, in a very 
simple setting. 

The paper is organized as follows. In section |I| we write the Fermion Boson model 
and the extended DMFT equations. Then we describe the semiclassical strategy for their 



analysis in both weak and strong electron-phonon coupling. In section |TT| we present results 
of solving the saddle point equations for 3d phonons coupled to electrons in different regimes 
and discuss the agreement with results in QMC approach. In section [IV] we describe the 
results for 2d phonons. If the electrons are fully integrated out the semiclassical treatment 
of EDMFT has to reduce to a mean field theory in classical statistical mechanics. In section 
[V] we compare EDMFT with other classical mean field treatments such as the Weiss mean 
field approach and the Bloch Langer method.@ The stability analysis of the EDMFT theory 



is carried out in Appendix [D]. 



II. MODEL AND SEMICLASSICAL APPROXIMATION 

The model under consideration is described by the lattice Hamiltonian: 

H = H e i + H ph + H e i^ ph (1) 

where 

H el = ~Yl tiAaCj* ( 2 ) 

HvH = £ || - E #•'•<•'•, (3) 

2 ij 



EAx^c^c^ - -) (4) 



icr 1 

The first term describes free electrons, c| CT (ci CT ) creates (annihilates) an electron with 
spin cr on a site i. The second term describes nonlocal (dispersive) phonons, Xj and Pi are 
canonical variables. The last term couples the fermionic and the bosonic degrees of freedom. 
We consider a half filled system of fermions. 

The second term could alternatively be written as: 

H ph = Y, UJ i( a l a i + h ( 5 ) 

q Z 

where a q , a q are related to the phonon field by x q = (2Mu q )~^ (a q + aL q ) and to^ = J q /M. 
Dispersion (momentum dependence) of the boson frequency uj q stems from nonlocal char- 
acter of Jij. The local limit of corresponds to the Holstein model0, so the model under 
consideration is an extension of the Holstein model to dispersive phonons. 

The extended DMFT equations for this mo dell are a set of equations for Weiss functions 
Go„(iuj n ) and D^^Un): 
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G a [G , D ](iu n ) = [ iuJ n -t q + G a [G , D ) 1 (iuJ n ) - G ^(iu n ) 1 (6) 

D[G ,D ){iu n ) = ]T [M{iu n f -J q + D[G ,D )-\iu n ) - D^iun)]' 1 (7) 
q 

where full Green's functions G a (icu n ) and D(iu n ) are expressed through GQ^(iu n ) and 
D^iiuJn) in terms of the effective impurity action: 

S eff = cl(iu n + iu m ) (G^{iu n )8o^ m + Xx(iu m j) c a (iuj n ) - -DQ 1 {iu m )x 2 (ko m ) (8) 



j^m iff 



n r \ JVH, Cq, x]c g (w n )4(w n )e s -n^M 

G a {iu n ) = — — _ s t - = (c a {iu n )ci{zu n )) Seff (9) 

D a (iuj n ) = -(x 2 (iuj n )) Seff (10) 

Eqs. ([^- |T0|) in general have to be solved numerically. In many cases though, a number 
of approximations reducing numerical work but preserving a physical content of the prob- 
lem is possible. One of the approximations is in using a model density of states (DOS) 
for fermions and bosons, so that momentum summations in EDMFT equations could be 
performed analytically. It is convenient to chose semicircular electron DOS: 

Pei(e) = ^- 2 VW^ (11) 

where W is the electron band halfwidth. The particular choice of semicircular electron DOS 
is qualitatively unimportant since we consider half filled electron band. For phonons, on 
the contrary, the shape of the phonon band near the bottom is crucial for temperatures 
smaller than the phonon band width. For d-dimensional phonons the bottom of the band 
has e^ 2 singularity. That is why to represent 3d and 2d phonons we chose semicircular 
and step-function-like phonon DOS respectively: 

3d p ph (t) = —2\ful - (e - ^o) 2 (12) 

71UJ v 



P P /^) = ^ H - (e - c ) 2 ) (13) 



2d , 

2cji 
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After replacing summations over wave vector by integrations over energy, Eq. (|^) and 
Eq. (0) read: 

G a (iuj n ) = J de^ (14) 



D(iu n ) = I dep^ (15) 



where ( = iu; n + (j~ 1 (iu; n ) — (ia> n ) , £ 2 = M(iu; n ) 2 + D _1 (zu;„) — DQ l {iuj n ); density of states 
p(e) = For electron p e i(t) and phonon p p h(t) DOS respectively e g = t g and e 2 = J g . For 
DOS defined in Eqs. (pTT| - p~3D integrations over energy in Eqs. ( pT4| , p~5D yield: 



G^ n ) = W2 [C-s^e-W^ (16) 

where s=sgn[Im£]. 

3d D(tcu n ) = -L ( 2 £ + ^-^o) 2 -^ - V^ + ^o) 2 -^ 2 ) (17) 

"(e + ^i) 2 -^, 



2d D{iuj n ) = — ^— In 



(18) 



.(e-^i) 2 -^ 

We consider here a semiclassical treatment of the problem. In its most general form, 
the approach has been described in ref,0 and is an application of the saddle point method. 
In this paper we use a more limited form of this method that consists of evaluating Eqs. 
(§,[K]) by a saddle point technique. It can be viewed as a combination of two separate 
approximations: the static approximation (equivalent to the phonon mass M — > oo limit) 
and a saddle point analysis of the EDMFT equations in the static approximation. 

The approach of treating the collective excitations as classical, while the electrons are 
treated fully quantum mechanically, goes back to the Hubbard approximation.0 It was 
pointed out that a static approximation of the impurity model coupled with the DMFT self 
consistency conditions indeed gives a solution closely related to Hubbard's.0 This approach 
has been used extensively in refsl^E in DMFT studies of the Holstein model. From the 
DMFT studies of the Mott transition,© we know that this approach becomes insufficient in 
the correlated metallic regime at very low temperatures, where a quasiparticle feature forms 



in addition to the spectral features produced in the semiclassical approximation. It is worth 
pointing out, that improvements of the static or of the saddle point approximation, will 
not remedy this shortcoming, which requires a non perturbative resummation of instanton 
events. Still, we show here that this simple analysis is able to reproduce all the trends of 
the solution of the EDMFT equations by the more expensive QMC method.@ 



The EDMFT equations in the static approximation Eqs. fl9|, |T0|) reduce to: 

G(iu n ) = I dxP(x) \ (19) 

J G Q (lU n ) + AX 

D = -pf dxP(x)x 2 (20) 



where 



P{x) = i exp L £ In (l - G^uj n ) 2 \ 2 x 2 ) - ^D, l x 2 ] (21) 



n>0 



Eqs. fll9| - pTl) have to be solved together with Eqs. (0, 0). In the static limit only the zero 
phonon frequency survives, so we drop frequency index for the phonon correlation functions 
Do and D . In Eqs. ( |19| - |21| ) and everywhere below we consider x being the phonon field 
amplitude, it is related to its Fourier transform as x = /5~2x aJm=0 . We consider no symmetry 
breaking in the electron spin channel, so we dropped the spin index; factor g (equal 2 for 
spin one-half) in the Eq. fl2~T|) appears from trace over the spin index. N normalizes P(x) 
to unity. P(x) is the probability distribution function of the phonon field amplitude x. 

We now evaluate Eqs. ( |T9"| - PTD in the saddle point approximation in the variable x. There 
are two limits, weak and strong coupling. In the weak coupling the saddle point is at x = 0, 
and in the strong coupling there are two equivalent saddle points at x — ±a;o 7^ 0. Deriving 
the saddle point equations we explicitly use semicircular electron DOS, Eqs. (|TT] , [161) . The 



relation between the bare and full Green's functions is especially simple in this case: 

GoiitUn)' 1 = iuj n - t 2 G(iuJ n ) (22) 
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where t = W/2. Everywhere below in the paper energy is measured in units of t. In this 
paper we restrict ourselves to the particle-hole symmetric case. 

In the weak coupling regime in the saddle point approximation, which includes Gaussian 



fluctuations of x around zero, semiclassical EDMFT equations Eqs. (|T9|-f21|) read: 



G(G + u) - (G + u) +a —0 (23) 

n> {G + u n y 

where G = iG(iu n ) and a 2 = X 2 \D\T = -X 2 T J dep ph (t)[D~ l - D Q X - e 2 ]" 1 , so a 2 is solved 
for the phonon self energy thus making the system of the saddle point equations closed. 

In the strong coupling regime we consider two saddle points x = ±Xq. We discard 
fluctuations around these points (so \D\ — (3xq), since nontrivial information is contained in 
the fact that we have two saddle points, and not in the Gaussian fluctuations, like it was in 



the case of weak coupling. EDMFT equations Eqs. ( pj| - |2l|) now read: 



G(G + oj) 2 - (G + u) + Go? = (25) 



D„ 1 = -T V —~ ^ (26) 

n Vo G + co n ) 2 + X 2 TD 



Weak coupling equations Eqs. (|23|j2^) are a saddle point expansion up to the first order 
in small parameter X 2 DT, and strong coupling equations Eqs. (|25| , |26"D - up to the first 



order in large parameter X 2 D/T. These equations have overlapped regions of applicability, 
provided Tel. This allows us to combine weak and strong coupling equations into a 
unique set of equations, controlled by the small parameter T: 

G(G + u) 2 — (G + u) + Ga 2 = (27) 

Do 1 " D- 1 = -2gX 2 T £ <?— (28) 

These are our final semiclassical EDMFT equations. They are exact in the limit MT 2 3> u)q, 
Tel. For 3d and 2d phonons Eqs. ( p7tp8l) have to be solved together with Eq. ( |l7l , |l"8"| ) . 



where £ 2 = D 1 — D . Saddle point equations Eqs. ( |27| , ^8D are very simple, they can be 



solved for D and D with minimal numerical efforts. Left hand side (lhs) part of Eq. 
is a third degree polynomial, so electron Green's function can be written as an elementary 
function determined by a single parameter a 2 which is a function of phonon self energy and 
bare parameters of the model. 

In the limits of small and large a 2 (or A) Eqs. flSTJZBD are solved completely for self 
energies: 

a < 1 



E d («*0= (-^ + \A + (^) 2 )« 2 (29) 



\h = -^A 2 (30) 



Moreover, in the dispersionless case a 2 = X 2 T/(uJq — ffA 2 ). This expression is valid every- 
where except for the small region A A ~ ujqT below A c ~ ujq. We consider here a disordered 
phase solution. In d = 3 the disorder solution becomes unstable at A ~ c^o — while it 
remains stable for all coupling in d — 2. The self energies in the strong coupling regime 
a ^> 1 are given by : 

a 2 

T, e i(iu n ) = — (31) 

= -€ < 32 > 

In the strong coupling the phonon field distribution function is split in two peaks. The 
peak separation is 2xo, where xq = —g\D /2, D = —{5x\. In the dispersionless case a 2 = 
(|A 2 /u; 2 ) 2 , this is valid when A ^> luq. This is completely similar to the previous analysis.0 

In d = 3 the instability to the ordered phase occurs already at small uj\ ~ j^x?, so Dq 
stays practically unrenormalized. 

In d = 2 at u\ ~ -j^zj? the system enters a regime when the phonon energy is expo- 
nentially small: 



Zph - (^o - ^i) 2 ~ 2o»i exp 



(33) 



In the limit T — > one readily obtains the polaron formation condition, which happens 
at intermediate (A c ~ u ) coupling: 

- ^D = 1 (34) 

where D = luq 2 in the dispersionless case, but has to be found numerically for interacting 
phonons. 

III. 3D PHONONS 

In this section we compare our semiclassical solution to the exact QMC results. The 
saddle point equations we derived are exact when (2-kT) 2 Muq 2 ^> 1 and Tt^ 1 <C 1. The 
QMC resultsi however, were obtained for (2ttT) 2 Muj 2 « 2.5 and Tr 1 « 0.13. We want to 
show, that even in these cases when the parameters controlling the saddle point equations 
are relatively close to 1, the semiclassical solution, even without including the refinements 
outlined in re 10 not only captures all the qualitative trends of the exact solution, but in 
many instances is quantitatively close to it. 

We study the case of 3-dimensional phonons. We use the same parameters as in the rei@: 
inverse temperature (3 = 8, the phonon band is centered at ujq = .5, electrons have double 
spin degeneracy g = 2 and hopping amplitude t = 1, phonon mass M = 1. The electron 



band is half filled. To model 3d phonons semicircular DOS Eq. (|12"D is used. We present 
the solution of Eqs. (0,H) and Eq. fll?]). 

In every figure in this section we plot both our and QMC curves. Our results are plotted 
using solid or dashed lines only, without symbols. QMC graphs are presented using dotted 
lines and always with symbols. 

A local instability, starting from the disordered phase, takes place within EDMFT when 



as discussed in refQ the effective phonon frequency to* = J (cu — uji) 2 + 11, given by the 
pole in the phonon Green's function, becomes equal to zero. The phonon mode softening 

9 



for different values of the phonon dispersion are shown in Fig.|l|. uo* is plotted versus the 
quantity characterizing the effective interaction: U = A 2 ju>\. The effective electron-electron 
interaction, mediated by phonons is given by (7 e // = A 2 _Do evaluated at zero frequency . 
When the phonon dispersion vanishes U e ff = U . The upper curve in Fig.[TJ corresponds to 
the dispersionless case. 

On the other hand, we find this is not the best way to detect an instability to an ordered 
phase and we discuss in Appendix [D] an alternative way to compute within EDMFT the 
phonon self energy which retains momentum dependence. 
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FIG. 1. d = 3. Effective phonon frequency uj* as a function of A 2 /w 2 at w i = 0.0,0.1,0.2,0.3. 
Comparison to QMC. 



The equations in section || do not include possibility of the phonon field symmetry 
breaking. They need to be modified to describe long range orderingji and we implement 
this in section 0. A well known property of mean field theories is that they allow the 
analytic continuation of mean field solutions beyond the parameter regime where they are 
stable. This was very fruitful in the understanding of the paramagnetic Mott insulating 
phase which is unstable to ferromagnetism.0 As was done by QMC in ref ,0 we study the 
continuation of the EDMFT equations beyond paramagnetic phase. It may hopefully be 
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understood as a metastable phase. This requires some care since the instability to a charge 
ordered phase is signaled by a singularity appearing in the integrand in Eq. (|TJ) and this 
instability causes D to acquire an imaginary part. As in ref,0 we take the principal part 
of the integrand, the imaginary part of Dq, being equal to zero in every numerical iteration 
loop, which allows us to compare our results with the results of the QMC method. 



A. Weak coupling 

The finite dispersion treated within DMFT renormalizes Dq (see Fig.Q(a)). Since the 
effective electron electron interaction is proportional to Do, the electron self energy is en- 
hanced as well (see Fig.^|(b)). While the features of the exact solution are qualitatively 
well reproduced in the semiclassical approach, it lacks quantitative agreement. The weak 
coupling is the worst case. The quantitative agreement is better for intermediate and strong 
coupling. 



11 




o 

1 2 3 4 5 6 



FIG. 2. d = 3. Weak coupling A = .2. Comparison to QMC. a) The bare phonon Green's 
function, b) The imaginary part of the electron self energy with the spectral function in the inset. 
uji = 0.1,0.2,0.3. 

B. Strong coupling 

In the strong coupling regime as dispersion increases, D renormalizes downward (see 
Fig.j3|(a)), together with the electron self energy (see FigJ^(b)). For the strong coupling the 
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quantitative agreement with QMC is very good. 
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FIG. 3. d = 3. Strong coupling. A = .8. a) Bare phonon Green's function, b) Imaginary part of 
the electron self energy with the spectral function in the inset. u\ = 0.1,0.2,0.3,0.4 

C. Intermediate coupling 

At intermediate coupling system is in a crossover between weak and strong coupling 
regimes. As increases, effective electron-electron interaction first becomes stronger, D 
and electron self energy increases, like at weak coupling. At uj* = behavior changes on 
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reverse, the picture is similar to strong coupling case. This is illustrated in Fig.[|(a) and 
Fig-l(b). 
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FIG. 4. d = 3. Intermediate coupling A = .4. Comparison to QMC. a) Bare phonon Green's 
function, b) Imaginary part of the electron self energy with the spectral function in the inset, 
wi = 0.0,0.2,0.4. 
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IV. 2D PHONONS 



In the previous section we calculated various functions at different parameters in the 
3d case. The saddle point approximation is exact in the limit of infinite mass M and 
zero temperature T. At finite M and T the applicability of the method in a wide region 
of parameters was established in the previous section by comparison to QMC data. In 
this section we study 2d phonon case (Eqs. (|2"7P%|) and Eq. flIS|)) in the same range of 
parameters. 

Unlike the 3d case, the 2d disorder solution is locally stable, and we focus on this solution 
in this section. As we will show in section |V| EDMF in dimensions d < 4 , gives rise to a 
first order transition at a critical coupling strength. We study the disordered state solution 
continued along the second order transition branch, which is skipped in the first order 
transition. 

In the EDMFT approach, d = 2, appears as a lower critical dimension for finite tem- 
perature second order transition. This result describes accurately the situation with order 
parameters posessing a continuous symmetry, but it is a spurious consequence of the inabil- 
ity of a local approximation to generate spatial anomalous dimensions in the cases where 
order breaks a discrete symmetry. 

First we illustrate the exponential softening of the collective mode : in Fig.[5] we plot 
effective frequency u* versus U . The Fig.[| should be compared to Fig{l| (3d case). In the 
latter the curves hit U axis, what implies second order transition. In Fig.[| the curves rather 
gradually approach U axis, never crossing it. 
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FIG. 5. d = 2. Effective phonon frequency lo* as a function of A 2 /ojq at 
wi = 0.0,0.05,0.1,0.15,0.2 

Phonons generate effective electron-electron interaction ~ A 2 Do , so we are especially 
interested in D behavior. We investigate 2d system for similar sets of parameters as we did 
in 3d case. 

We obtained the following plots in weak, intermediate and strong coupling regimes: 
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FIG. 7. d = 2. Intermediate coupling. A = .4. a) The bare phonon Green's func- 
tion, b) The imaginary part of the electron self energy with the spectral function in the inset, 
wi = 0.0,0.1,0.2,0.3,0.4 
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FIG. 8. d = 2. Strong coupling. A = .8. a) The bare phonon Green's function, b) The imaginary- 
part of the electron self energy with the spectral function in the inset. u>\ = 0.0,0.1,0.2,0.3,0.4 
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FIG. 9. d = 2. Bare phonon Green's function. A = .2, .3, .4, .8 

The behavior of 2d system is very similar to that of 3d system before the energy of the 
phonon mode vanishes. In all cases Dq gets renormalized, as the dispersion and consequently 
the effective interaction increases. The electron self energy enhances correspondingly. The 
only difference is the rate of D renormalization in the weak, strong and intermediate cou- 
pling. D renormalizes faster at larger A (see Fig.|3|), since electrons are stronger coupled to 
phonons. 



V. ORDERED PHASE AND CRITICAL TEMPERATURE 

We now turn to the generalization of the EDMFT equations to the ordered phase.i For 
simplicity we will consider a classical model. This is justified, since in the semiclassical limit 
we can always integrate out the electrons, reducing EDMFT equations to classical mean field 
equations. For instance tracing out the electrons and performing a static approximation in 
the electron-phonon field leaves us with an action of the form ( neglecting terms of order 
(f) 6 and higher): 
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sw = E \ti + jti - E *y = SloM] " E 

t ij ij 



(35) 



To extend the EDMFT approach to the ordered phase it is useful to write down the 
Baym Kadanoff functional for the action, 



T[m,D] = --Tr\ogD + -TiD Q 1 D + -mD 1 m + <S>lm,D} 
2 2 2 



(36) 



$ is a sum of all two particle irreducible diagrams constructed from phonon Green's functions 
D, phonon field expectation value m and four legged interaction vertex —3\U. We could also 
say that $ is a sum of all two particle irreducible diagrams constructed from phonon Green's 
functions D and four, three, two legged vertices plus the first diagram shown in figure (|IT)|). 
which contains no propagators. The vertices yield factors of —3\U, —3\Um, —3Um 2 and 
— jUm 4 for four, three, two and zero legged vertices respectively. Each diagram in $ has an 
extra —1 factor. 



In the Fig.|10]we drew first and second order (in U) diagrams entering $. The extended 
DMFT equations in the ordered phase are derived by making the local approximation on $ 
in the Baym Kadanoff functional and solving the stationary conditions for the magnetization 
and the local propagator resulting from the stationarity of Eq. fl3"E| ) after this local approx- 
imation is made. In the local approximation the leading terms in a perturbative expansion 
in the quartic coupling are given by $ = \Um 4 + \XJDm 2 + | UD 2 — 3U 2 D 3 m 2 — jU 2 D 4 + ... 



+ 




+ 




+ 



FIG. 10. Diagrammatic expansion of $ up to the first two orders in U 



Stationarity of the functional in Eq. (|36D would give exact equations for D and m. In 
the local approximation these equations reduce to EDMFT equations in zero magnetization 
and therefore generalize those to the ordered phase.! They are given by 



m(r - J q=0 ) + Sm 
21 



D = E[r~J q + 2^ (37) 

where only local graphs are included in $ . Diagram series equivalent to the first equation 
are shown in Fig.|ll|. 




• = h» + *- + h» + 




FIG. 11. Diagrammatic expansion of m up to the first two orders in U. A thin line is the free 
phonon propagator Dq, a thick line is the full phonon propagator. A full dot stands for m. 

For the practical solution of the EDMFT Eqs. ([37]) , it is useful to follow the dynamical 
mean field procedure of introducing an impurity local effective actionEUil to sum up the 
graphs generated by the functional for $ and its functional derivatives and 5^ /8m 

in terms of the cavity fields h and A . The effective action 

Sedmft{4>] = Si oc {4>] -h(j>- y0 2 (38) 

generates the correct local quantities provided that the Weiss fields h and A are chosen to 
obey the EDMFT self consistency conditions: 

r — A — D = -2°— (39) 

h = -m(r - A) - — (40) 
dm 

Eqs. ( 1371) and Eqs. (|39| , |40D are a closed set of EDMFT equations, describing both ordered 



and disordered phase of the classical system Eq. (|35|) : 

D = J2[D- 1 + A-J q ]- 1 (41) 



Sedmft[4>\ = Si oc [4>\ - m (J q=0 - A) - y0 2 (42) 
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m = ((p) s 



(43) 



D=(^)s EDMFT -m 2 (44) 

The equations above are consistent in describing the transition: magnetization vanishing 
in the ordered phase, and divergence of spin susceptibility across the transition do occur at 
the same critical temperature. 

When U — > oo, r — > — oo, U/r = — 1 the system, described by the action Eq. (|35"D 
reduces to a classical Ising model with spin values ±1. In this limit the standard Weiss 
mean field equations 

m = tanhmJ g=0 (45) 

can be compared with the EDMFT equations which now read: 

m = tanhm (J q= o — A) 
1 - m 2 = £ [(1 - m 2 y l + A - J q y l (46) 

Q 

We will also compare the EDMFT equations to an extension of mean field theory due to 
Bloch and Langer (BL)J3 

M x = f dx{2nG 2 )-^ eX p[- (Ml Jq= ° - % ^ } tanhx (47) 

J 2&o 



M 2 = [ dx(2nG 2 )-^ exp[-^%^— ^] (cosh x)~ 2 (48) 
J 2G 2 

— oo 

G -?T^ (49) 

where Mi is the magnetization. It can be shown that EDMFT counts (without overcount- 
ing) more terms in diagrammatic expansion of various physical quantities, like correlation 
function or free energy, than BL method does. One can also check that for the classical ±1 
spin model EDMFT gives a better estimate for T c than BL method does. 
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We computed the critical temperature T c for the Ising model on a Bethe lattice with 
finite coordination following the paramagnetic solution till it disappears using the different 
approximations described in this section. The results are shown in Fig.[L^. The EDMFT re- 
sult shows significant improvement over MFT, and is slightly better than BL method. Some 
more technical details comparing the approximation schemes are relegated to Appendix |A]. 
0.6 




FIG. 12. Top to bottom: BL, EDMFT, exact solution, MFT; Critical J c vs nearest neighbors 
number z 



In spite of the quantitative improvement of T c , the order of the transition is incorrectly 



given by the EDMF approximation. In Fig.[13] we present result of solving EDMFT equa- 
tions for the simplest ±1 spin model in d = 3. We plotted magnetization as a function of 
temperature. At sufficiently low temperature the solution consists of three branches with 
magnetizations m 3 > m 2 > mi = 0. These branches are extrema of the free energy, which 



is shown schematically in Fig.|14]. m 3 and mi correspond to local minima in the free energy, 
while 1712 corresponds to a local maximum and is unphysical. The transition is clearly of the 
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first order. The order of the transition does not change up to d = 4. 

The inability of EDMFT to predict the correct order of the transition is related to 
the inability of a local theory to produce anomalous dimensions, and persists in quantum 
problems when the dynamical critical exponent and the dimensionality are such that they 
require the introduction of spatial anomalous dimensions. Details are given in Appendix [B| 
and 0. 



1 



m 












FIG. 13. Magnetization vs. temperature on a cubic lattice. There are three branches at 
3.96 < T < 4.45: 777,3 > rri2 > m\ = 0, m\ and r/13 are physical solutions, while 777-2 is not. 
Classical ±1 spin Ising model. 
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-1 1 

m 

FIG. 14. Free energy evolution with T. Free energy has a single minimum (m = 0) above 
T = 4.45; at T = 4.45 the solution bifurcates at m = .45. There are free extrema at 3.96 < T < 4.45 
corresponding to = ml < m2 < m3. As T aproches T c = 3.96, ml merges with ml. 

VI. CONCLUSION 

We have performed a semiclassical analysis of the EDMFT equations for a simple Fermion 
Boson model. Comparison with earlier QMC treatments of the same problem reveals that 
this method reproduces semiquantitatively all the trends found in the previous study! It 
can be used therefore in the study of more complicated systems, such as fermions interact- 
ing with spin fluctuations. We have also investigated this approach in the ordered phase 
revealing some inadequacies of the approach which are closely related to the existence of 
anomalous dimensions in finite dimensional systems. Since this non trivial k dependence 
which is characteristic of low dimensional systems cannot be generated by a local theory, 
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EDMFT produces spurious results such as the existence of a first order phase transition in 
d < 4 . Since at zero temperature the dynamical critical exponent is such that in two dimen- 
sions an ansatz without anomalous dimensions is internally consistent ,0 a continuation of 
the disordered state, beyond the first order phase transition, might be useful to study this 
system. In this spirit we pointed out that a continuation of the EDMFT, at finite tempera- 
tures, suitably interpreted, gives improved estimates of the critical temperature compared to 
the simplest mean field treatment or the Bloch Langer method.0 It could be used to obtain 
better estimates of ferromagnetic transition temperatures where spatial fluctuations of the 
order parameter substantially decrease the Curie temperature below the DMFT estimates. 
This is the case of bcc Iron,@ a problem which will require a more realistic investigation 
of EDMFT. Further investigation of the quantum problem will require zero temperature 
methods which go beyond the semiclassical approximation. 



APPENDIX A: CRITICAL TEMPERATURE 

The transition to the ordered state for a classical model in BL method is signaled by the 
divergence in the zero momentum term in the sum: 

G 2 = y — ^ — (ai) 

This equations is analogous to the selfconsistency equation of the EDMFT, it is arises when 
summing ring diagrams. M 2 is a vertex and Jjj is a line in the ring diagram. G 2 (i)M 2 is the 
sum of all ring diagrams which cover the site i. M 2 and G 2 are related to D in EDMFT: 

D = M 2 + M 2 G 2 M 2 (A2) 

Below we are explicitly summing ring graphs on a Bethe lattice to express G 2 through 
M 2 . 

We are introducing notations, G 2 (i) and Q(i): G 2 (i) equals G 2 (i) when the latter is 
computed on a lattice with all but one bonds cut out from the site i. Q{i)M 2 includes those 
diagrams from G 2 (i)M 2 which have only one vertex belonging to the site i. 
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The following relations can be established: 

G 2 = zG 2 + zG 2 (z - l)M 2 G 2 + zG 2 (z - l)M 2 G 2 {z - l)M 2 G 2 + 



(A3) 



G 2 = Q + QM 2 Q + QM 2 QM 2 Q + ... (A4) 

Q = JM 2 J + JM 2 (z - l)G 2 M 2 J + JM 2 (z - l)G 2 M 2 (z - 2)G 2 M 2 J 

+JM 2 (z - l)G 2 M 2 (z - 2)G 2 M 2 (z - 2)G 2 M 2 J + ... (A5) 

where z is the number of nearest neighbors, J is a bond on the lattice. 
Summing geometrical series, we obtain: 

G 2 = ^ (A6) 

1 - (z - 2)M 2 G 2 

G 2 = ^ (A7) 

1-M 2 Q K } 



Q = .PM 2 1+M ' G '_ (A8) 



Solving these equations we get: 

(l " ~~ 1 - zM 2 Q 



G 2 = - %— (A9) 



(z - 1)2M 2 
Eq. (|£|) and Eq. ( |AT0l ) solve G 2 for M 2 



1- a/1- (z- 1)(2M 2 J) 2 

q = — — ~ — — ~ — — (aio) 



Curves J c vs z for BL and EDMFT, together with MFT solution J c — 1/z are presented 
in Fig.[[| and compared to the exact solution: 
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APPENDIX B: ORDER OF THE PHASE TRANSITION 



Here we are proving that the EDMFT equations give a transition of the first order for 
d < 4 and of the second order in the higher dimensions. For classical phonons EDMFT 
equations read: 



mir — J,, 



>q=0) 



£ = E[r + 2 «M_ jj-i ( Bi) 
i 



It is easily seen that 



5 2 $[m,D] 



5m 2 



= 2^M1 (B2) 
5D y 1 

m=0 



Solving Eq. ( pT|) for r using the above relation for derivatives, up to the second order 
in m we have: 



2 

m 



-^--D |m=0 — D\m=0 — ^2 



(B3) 



The coefficient in front of m 2 in the right hand side (rhs) is positive. The lhs of Eq. ( B5[) 



is oc m 2 ; while the rhs has two contributions, one oc m d ~ 2 and the other cx 5(3, where 
5(3 = f3 — (3 C . For d < 4 the term oc m d_2 is dominant and 5(3 oc —m d ~ 2 < 0. A negative 5(3 
implies the first order transition. For d > 4 the term m 2 from lhs becomes dominant and 
5(3 oc m 2 > 0. This is the usual mean field behavior resulting in a second order transition. 

We showed that in a classical model the transition is of the first order below the upper 
critical dimension. The same is true for a quantum transition as well. We show it in the 
appendix considering large iV limit. 

As discussed earlier in connection with the order of the transition, this artifact of the 
EDMFT results from the inability of a local theory to capture physics that requires the intro- 
duction of anomalous dimensions. In spite of this shortcoming, when properly interpreted, 
EDMFT results in improved estimates of the critical temperature relative to DMFT. 



29 



APPENDIX C: QUANTUM PHASE TRANSITION 



In this appendix we investigate the phase transition in the quantum version of 4 model. 
We compare EDMFT and a full lattice model using large N technique. We will show that 
above the upper critical dimension d > d uc = 4 — z the exact critical exponents and the 
critical exponents obtained in EDMFT coincide. Below d uc the EDMFT and the lattice 
model exhibit different critical exponents. In the EDMFT the transition is of the first order 
for \d uc < d < 2 and of the second order otherwise. The transition is of the second order in 
the lattice case. Moreover, in EDMFT the exponents have universal value for d < \d uc and 
a non universal value for \d uc < d < d uc . 

The lattice model is described by the action: 

s = \d?? + *Ltff (ci) 

where D^} q = r + \u\ * + g 2 , 2 = J2a=i 0a; U = u/N, r is & variable parameter which drives 
the phase transition. Corresponding EDMFT equations are: 



'Ou,q=0 



+ ' ' J = (C2) 
dm 



»=K + ^1- (C3) 

The functional $[m, D] includes all two particle irreducible diagrams which are con- 
structed from the magnetization m (dot), the particle propagator D (line) and the interac- 
tion term U (four legged vertex). $ satisfies the following equation: 



5 2 <5>[m,D} 



<J$[0,£>] 

= z 

m=0 



2^ (04) 



Sm 2 

Expanding $ in small m and using Eq. flCjj ) we write EDMFT equations as: 

6*[0,D] U**[m,D] 2 _ 
D ^=° + 2 ~5D~ + 6 6m' lm ^° m "° (C5) 

D a = £dq d [\u\l +q 2 + {2T - l ^^ Mm 2 ]- 1 (C6) 
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where T = 

Let D 0c , D c , r c be values of D , D, r in the transition point. Subtracting 
-^ocL,g=o + 2^jjy^ = from Eq. ( |C5|) and keeping lowest order terms we have: 

5r + 2T5D + \^m 2 = (C7) 

where 5r = r — r c , 3D = D — D c . This equation provides a relation between the variation 
of the driving term r and the order parameter m. We will show that for d > d uc the last 
term in the left hand side wins over the second term, the transition is mean field like. The 
second term becomes important and determines the character of the transition for d < d uc . 

r = x + 

FIG. 15. 1/N expansion of T. All diagrams are of the order 1/N 

We will consider the large N limit up to the order of 1/N. Diagrams which enter F are 
chains of bubbles (see Fig.|T5|), which can be summed up as geometrical series: 

where x w ~ / di>D v D u+UJ in the quantum problem or x — D 2 in the classical problem. The 
only term of the order 1/N which enters £ 4 <I>/£ 4 m is 6u/N. All other terms are of order 
0(1/N 3 ). 

Eq. (1C7D and Eq. ( |C8D holds in case of a lattice as well, but D now depends on both 
momentum and frequency, and summations now run over wave vectors as well. 

The upper and lower critical dimensions are determined by the convergence of integrals 
TrSD and TrD in the ultraviolet and infrared limits respectively: 

Tr5D ~ f dudq % (C9) 

TrD ~ f dudq % (CIO) 

J (\uj\z+q 2 ) 

These equations are the same for the mean field and lattice models, they yield the upper 
critical dimension d uc = 4 — z and the lower critical dimension di c = 2 — z = d uc — 2. 
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We first consider EDMFT. In a crude way one can estimate: 

d > 2, D U ~J d d q{\u\i + q 2 )- 1 ~ (d - 2)~ l (A^ - \cu\^) (Cll) 

d<2, ~ —(d — 2) -1 |a;|^ 2 (C12) 

and 

d>2, X(tf ~JdvD u D u+u} ~ (d-2)- 2 H^ 2+1 (C13) 

d<2, x ^-(d-2)- 2 (2^ + l)- 1 |^| 2 ^ +1 (C14) 

z 

A g is a momentum cutoff. We see from Eq. (|C14j ) that for d < c? uc /2 the susceptibility \u> is 
divergent at low frequency, it leads to a universal critical behavior for d < d uc /2, as follows 
from the self energy calculation below. The self energy in large N limit is <5£ ~ 2TD: 

^<d<d uc , ^ ~ L(d - 2)- 1 u\l;\^+ 1 (C15) 
d lc <d<^, Z^±Jdvx- 1 D u+ ^l(d-2)(2^-l + l)\Lur^ (C16) 

In a similar way we can calculate a contribution from m to T5D 

^<d< d^, T5D ~ -^(d - 2)~ W" 2+ * (C17) 

d lc <d<^, T5D~}-td-2){2— + l)m- d+2 (C18) 
2 TV 2 



This result together with Eq. (|C7|) suggests that the transition is the first order for 



|rf uc < d < 2. 

Now we consider the lattice model. 

X., q ~ / d^ d p(k + ^l 1 + (P + ^T'dH 1 + P 2 Y l ~{d + z- ±y\\uj\i + g 2 )^" 2 (C19) 

If 1 2 

J dvdpx-lD u+u)<p+q ~ —(d + z-4)(\tj\-\n\uj\ + q 2 \nq) (C20) 

In this case the frequency dependent part of the self energy can be conveniently exponenti- 
ated to yield: D ~ Wuo^+q 2 -"!}- 1 with z = 2-N- 1 (d-d uc )c 1 (d) andr? = -N~ 1 (d-d uc )c 2 (d), 
where Ci(d) and 02(d) are some smooth functions of d. 

We also calculate a contribution from m to T5D ~ N~ 1 (d — d uc )m 2 In m. It yields 
<5r ~ with /3 = \ + (c? — d uc )j^C3(d). The transition is the second order in this case. 
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APPENDIX D: INSTABILITY ANALYSIS 



Let us consider a very general electron phonon Hamiltonian which describes an electron 
phonon system with electron-electron interaction (local or long range), electron-phonon in- 
teraction and phonon-phonon interaction (phonon unharmonicity) . We can always use a 
Hubbard Stratonovich decoupling on electron-electron interaction, so we assume that infor- 
mation about long range electron-electron interaction is stored in the phonon dispersion and 
we will not write the long range interaction explicitly. We can introduce a source dependent 
action S where the sources are coupled to different fields. The free energy W = — In / e s is 
the generating functional for expectation values of those fields. 



/ 



dxdx' 

cUx)G 0a \x - x')c a (x') + ^(x)D \x - x')<f>(x') 

+5(x - x ) {Un^{x)ni{x) + V A (f) A (x) + A0(x)4(x)c CT (a;)) 

-J„(x, x')cl{x)c a {x') - \<j)(x)K(x, x')(f)(x') - 8(x - x')L(x)(f)(x) (Dl) 

x variable includes both space and time in the above formula and repeated indices imply 
summation. Expectation values of the fields coupled to the sources are given by: 

G = 77- K = 2 W m = lL' (D2) 

Exact Green's functions correspond to the limit of zero sources. To study phase tran- 
sitions, like the transition when the phonon field acquires non zero expectation value, one 
needs to have the free energy as a functional of correlation functions only. Such functional 
can be derived as a Legendre transform of the free energy: T = W — JG — K/2D — Lm. 
The sources J,K and L have to be solved for G,D and m. The functional Y is called a 
Baym Kadanoff functional and its stationarity yields equations for zero source correlation 
functions. We present the functional without derivation: 

T BK [G,D,m] = 

TrlogG - Ti^ 1 - G- r )G - ^-TrlogD + l-TrD^D + ^-mD^m + <5>[G, D,m] (D3) 

Zj Zj Zj 
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G and D are free fields of the action, $ functional is the sum of all two particle 
irreducible graphs constructed from the original bare interaction vertices, from vertices gen- 
erated by shifting the phonon field by m and from full correlation functions G and D. 

The charge ordering instability can be studied by looking at the zero frequency- 
momentum phonon propagator behavior: the propagator diverges in a charge density wave 
(CDW) transition. Alternatively one can study the transition from the ordered side, by 
observing the order parameter vanishing (m in our case). The two approaches should give 
consistent results. We will first show that this is indeed the case in the exact theory, then 
we explain how a similar approach can be applied in the EDMF theory. 

Let us introduce some compact notations we are going to use. o a a is a field operator of a 
kind at a space-time point, a = G specifies an electron field operator and a = D specifies 
a phonon field operator. 0^ 1& is a four point function, which is a subset of all connected 
diagrams in the perturbative expansion of (o^ OpOg cfy , rules for selecting the subset of 
diagrams depends on a particular operator. Multiplication of two operators is defined by: 
[0^0^] a /3^s — Hfu> O^p ,ixv0^u>,-y5- We are introducing three four point operators: 1) Xo 
includes all graphs which enter skeleton graphs without interaction vertices. 2) £ includes 
all ID irreducible diagrams; 3) T includes all 2P irreducible diagrams. In our case reducibility 
of O a p^s is understood as disconnecting af3 part from ^5 part. "ID irreducible" means "one 
particle irreducible with respect to cutting a phonon line". "2P irreducible" means "two 
particle irreducible". Xo is trivially expressed in terms of correlation functions: x GD — 
Xq G = 0, X GG p n s — G ai Gp$ and Xoap^s = D ai Dp$ + D aS Dp 7 . 

We can write out the following Dyson equations for the components of E operator: 

S GG = X GG + ^GG^GG + ^GD^DG 
Yi GD = x GG V G GTi GD + X GG T GD^DD 

s DG = xo D r DD z DG + Xo DD r DG s GG 

s DD = xo D + x^r DD £ DD + x ? D r DG s GD (D4) 

or we could simply write: 
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(D5) 



Solving for £ we find: £ = [% 1 — T] 1 = — (<9 2 r BA -) x . Second derivative d 2 Y BK is 2 x 2 
matrix defined by: 



(d 2 T BK ) ab 



d 2 T 



BK 



(D6) 



OC a dC b 

where C is a two component vector: Cg = G, C B = D. 

£ matrix is related to the phonon self energy H p h in a simple way, as can be seen from 
the diagrammatic series in Figjl6|. 
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FIG. 16. Diagrammatic expansions for S and 

£ comprises all four legged 1PI graphs, while Sp^ comprises all two legged 1PI graphs. 
2PI four legged block is nothing but T. Two horizontal lines represent a couple of correlation 
functions of the same kind, GG or DD, we assume that a summation runs over each couple 
of horizontal lines, while 2PI four legged blocks are understood as 2 x 2 matrices. The first 
line is a diagrammatic analog of Eq. ( P5|) . The second line provides connection between E 
and Il p h- That can be written as: 



or in a slightly different way: 



+ 



to6 



dmdm dmdC, 



dChdm 



(D7) 



d 2 $ , d 2 T BK f d 2 T B K \ 1 d 2 T 



BK 



(D8) 



dmdm dmdC \dCdC ) dC'dm 

The condition for the CDW instability at wave vector q is: -D^ 1 — = 0. 
We will reproduce the above result studying CDW transition from the ordered phase. 
T B k is the free energy, so in the transition point 
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d 2 T 



IJK (D9) 



dmdm 

From the way Ybk is constructed it follows dTsK/dC = and dTsK/dm = dTsK/dm. If 
we use d(dT bk / dC) / dm = and Eq. ( P9|) we find: 

d 2 T BK d 2 T bk ( d 2 Y BK \ ^ d 2 T BK _ Q (mQ) 
dmdm dmdC \dCdC ) dC'dm 



This equation is identical to Eq. ( |L)8| ) as should be in exact theory. In EDMFT approach 
we take the local approximation for the two particle irreducible graphs. All 2PI graphs in 
r bk are contained by $. So the condition for m vanishing is still given by Eq. (|D10|) with 
$ being local. Alternatively we can use Eqs. (|D5|) where T ab = 9 ^ 2 g Cb is local, in which case 
these two methods are equivalent. Let us consider the second method, when the transition 
is approached from the disordered phase. 

The local r afe can be computed using the impurity action of EDMF theory. For simplicity 
we consider electron phonon interaction only, with the coupling A = 1. Equations similar 
to Eqs. ( p4| ) can be written for the susceptibility x ab = (o a ^o a o b ^o b ). In short notations it 
reads: 

X = Xo + Xof X (DH) 



where T is different from T of Eq. ( |D5| ) , because now it includes ID reducible diagrams. 
The relation between T and T is simple: 

f = r + A) (Di2) 

where D is 2 x 2 matrix, t)^ G = D and D^ G = D GD = D^ D = 0. Using Eq. (|E>D, Eq. 
( pil ) and Eq. (|D 1 2|) we can express the self energy S through the quantities which are 



directly computed from the impurity action: 

£ = \[XimpY l - [XoimpY 1 + D 0imp + [xo] -1 ] -1 (D13) 

where xtnp = {o\o a o\o b )i mp , x<K mp = &abC 2 a and D 0imp is the Weiss field of the impurity 
action: 
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S imp = J drdr' (D14) 
cl(r)G 0l l p ^r - t'K(t') + ^{r)D^ mp {r - r')0(r') + 5(r - t')0(t)4(t)c ct (t) 

The described method is exact in the limit d — > oo. At finite g? it yields a higher T c than 
a naive local approximation T, p h = 5&/5D. Assuming Y, p h = 5&/5D would be equivalent 
to taking V bk as being local in Eq. ( |D8| ) , while the correct approach is to take a local 
approximation on the $ functional only, not on the whole Baym Kadanoff functional. 
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